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ABSTRACT; 

Tne temporal and spatial variations of the index of refraction cause fluctuations 
of sound phase and amplitude that can be completely understood only by defining 
the index in terms of the duration, location, range and time of the acoustic 
experiment. A truncated ’’universal" spatial correlation function of the index has 
been derived from a simplified form of the Kolmogorov-Batchelor spectrum of 
temperature fluctuations in a homogeneous, isotropic medium. Although this 
correlation function is shown to be predictable simply from the depth of the 
experiment, it is of only limited validity with respect to large spatial lags. 
However, a Gaussian extrapolation of the "universal" correlation function 
together with the standard deviation of the index provides simple useful predic- 
tions of the sound fluctuations due to temperature microstructure in the upper 
ocean. 
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1 . Introduction 



Sources of Sound Phase and Amplitude Fluctuations 



As it proceeds away from its source, a sound at sea encounters changing 
values of temperature, salinity and density, as well as entrained objects and 
bubbles in motion. These acoustically significant quantities are constant 
neither in space nor in time. Slower than the microscale fluctuations due 
to surface waves and turbulence of the medium, which show perceptible changes 
during a time scale of the order of seconds , are changes due to internal waves 
with periodicities of the order of minutes , tidal and diurnal variations , and 
seasonal changes. Depending on the duration of the study, some of these 
fluctuations may appear to be simply periodic, some show the spectral char- 
acteristics of a narrow band noise, many are describable only by statistical 
methods , and still others , such as the appearance of schools of fish or the 
incidence of storms, are largely intermittent in character. (Weston 1969). 

It is important to make the point that it is not the inhomogeneities but rather 
it is the temporal change in the position or character of the inhomogeneities 
that is the source of interest in this report. 

The phase of a sound wave travelling through a stationary medium depends 
on the speed of sound along the path. This dependence is expressed by 
equations 1.1 and 1.2: 



velocity along the ray path 

k = wave number along ray path 
T = temperature 
S = salinity 
H = depth 

n(a)da = number of bubbles of radius, a, in increment, da 



The change of amplitude of a non-divergent , unattenuated sound beam 
along a fixed path through the water is caused by the lens-like patches or 
layers in the medium which create convergences or divergences of sound power. 
This focusing and defocusing changes with movement of the patches, and the 
sound amplitude at the receiver thereby fluctuates. 




(/5= J k ds = J 
<p= sound phase 



( 1 . 1 ) 



where c=c (T,S ,H» ft(a)da) 



( 1 . 2 ) 



|v| = magnitude of component of material 
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1 . 1 Microstructure in Mixed Water 



It will be many years before the turbulent, patchy, nature of the ocean 
will be unscrambled by fluid dynamicists and physical oceanographers. How- 
ever, a few general concepts and landmark ideas of turbulence theory will be 
helpful in guiding our thinking particularly about the upper ocean. 

Near the ocean surface, depending on the time of day and the meteoro- 
logical history at the location, there is often a well-mixed layer of so-called 
"isothermal" water. The name is a reflection of the fact that common bathyther- 
mographs (BTs) , show such a region as isothermal. In fact, temperature probes 
of higher resolution and shorter time constants have revealed that the crudely 
isothermal region is characterized by drifting patches of water of slightly dif- 
fering temperature. Fig. 1. 




Fig. 1 Temperature fluctuations Run 6,4.3 Meters, 0354, 22 October 1971 

Within Such a region there are temperature fluctuations at a point whose 
standard deviation, o[^, may be of the order of 0.1 C. We will assume that 
this value of cfj, may be obtained at a point by averaging over an appropriate 
time, or at a given time by averaging over an appropriate local space 
(ergodicity) . These values are a function of time of day and depth of measure- 
ment. The short-time standard deviation of the temperature fluctuations meas- 
ured with respect to the mean temperature is generally greatest near the surface, 
greatest during the late afternoon, particularly during sunny periods. (Sagar 1960) 

For example, measurements (Seymour 1971) made during the late afternoon 
of a sunny day in October, one mile off San Diego, Calif, reveal a decrease 
in otj, from 0.15to0.14to0.13C as the measurement depth was increased 
from 7 to 9 to 14 meters. The effect of night time cooling is interesting; by 
early morning of the following day the rms temperature fluctuations in the 
experiment had decreased to approximately 0.06C° and showed no clear 
dependence on depth. The water temperature gradient was only about 0. 1 0/ 
meter, being warmer at the surface. 

Although material velocity variations generally have a lesser acoustic 
effect than temperature variations (Neubert 1970), there are circumstances 



5 



where the motion of the medium is significant, and may even overwhelm the 
temperature influence. Two situations should be considered: Near the surface 
of the ocean the underwater material velocity will be almost totally dependent 
on the surface wave height. For long-crested (swell) waves in deep water 
the vertical and horizontal, waveward, components of the material velocity for 
any frequency component of the surface wave will be almost totally orbital 
"near" the surface. 

The depth dependence is given by 

V (x^,Z) = V (x^) e ^ w^ (1.3) 

where x = 2T1/A = surface wave propagation constant or 

w w , 

wave number 

Z = de pth 

V (x ) = orbital velocity of the x wave component at the 
w , 

surface 



Let us call the above underwater motion that is completely, and directly, 
ascribable to the surface wave action, the "coherent component" of the 
material velocity. Superimposed on the coherent motion of the medium there 
will be the "incoherent" turbulent motion, often, and for simplicity, assumed 
to be independent of direction. The region in which three dimensionally 
isotropic, or two dimensional horizontally isotropic, turbulence can be assumed 
to exist is a function of depth and of sea surface wave conditions. As a 
second example there is evidence (Dunn 1965) of horizontally isotropic turbu- 
lence in Menai Straits, 30 minutes before high tide, during low sea states at 
depths of 10 feet and below. The material velocity was not isotropic at 5 ft. 
de pth . 

Since the temperature fluctuations in the microstructure at sea are generally 
very small compared to the average temperature it is appropriate to use the 
linearized form of the speed dependence on temperature. 



c = c + bA T 
r 

Ac = b(T) AT 

where c^= reference speed at temperature, T (1.4) 



Ac 

b(T) 



c 



c 

r 



AT 



linear rate of change of speed with 
temperature at T 



(1.5) 
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1 . 2 Layered Micro structure 



It is now recognized (Woods 1966) that there are large regions of the 
ocean within and below the thermocline , which, for extended periods of time, 
possess an ordered structure of alternate layers of turbulent and relatively 
stable water. About one half of the ocean's regions below the mixed upper 
layer may show this very extensive active layered character. The active regions 
contain identifiable "layers", warmer than the liquid above or below them, 
ranging from thickness sometimes as little as 10 cm or less and extending for 
hundreds of meters, up to layers of thickness 10 to 15 meters extending for 
tens of kilometers. The layers are separated by thin "sheets" of relatively 
steeper temperature gradients of order 0.01 C /cm. Temperature inversions 
of the same steepness with compensating salinity variations are found in these 
active regions . And scattered throughout the active volume there are intermittent 
and patchy regions, again, larger in extent than in thickness. It appears that 
(Nasmyth 1970) turbulence microstructure in these layers is always accompanied 
by temperature microstructure . On the other hand temperature microstructure 
has been found that it is virtually free of turbulence microstructure , and that 
may represent evidence of a former turbulent microstructure. 

1 . 3 Other Underwater Structure Features 



Frontal regions of moving water masses have been found, paralleling the 
meteorologist's well-developed picture of air mass behavior. For example, one 
might expect (Nasmyth 1970) a front with an angle of advance of 1.5 , in which 
quiescent warm water is being uplifted by a mass of highly turbulent colder 
lower water. It is not yet clear how common such situations are in the oceans 
of the world. 

Internal waves, generated within the ocean mass, are the analogue of 
ocean surface waves. They are volume gravity waves which possess their 
maximum vertical displacement amplitude at the interface between two water 
masses of different densities, but their motion is detectable far above and 
below this interface. All the characteristics of the water mass will change at 
a point through which internal waves pass. For example, temperature fluctua- 
tions of a few degrees with periodicities of the order of minutes or hours are 
commonly found. 



2. Spectrum and Spatial Correlation of the 
Index of Refraction 



Studies of temperature and velocity structure within the ocean concentrate 
on the wave number spectral description or its transform, the spatial correlation. 
It is clear that patchiness of salt concentrations or bubble content at sea could 
also be described effectively in these same terms. The goal in turbulence 
research has been to find a universal spectrum. 
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2 . 1 The Wave Number Spectrum 



The underwater material velocity that is uncorrelated with the ocean 
surface displacement, the turbulent velocity, cannot be completely isotropic. 

It is easy to see that the very long wave lengths of the Z component, at least, 
would be affected by the depth of the water and the depth of the mixed layer. 
Nevertheless , it is convenient for the theorist to talk of localized isotropic 
turbulence as an idealized situation. 

The guiding light in turbulence studies has been the Kolmogorov 
Hypothesis which deals with isotropic turbulence. The Hypothesis states that, 
regardless of the mean flow and large scale motions (long wave length) from 
which turbulent energy comes, the small scale components of a "fully 
developed" turbulent region are in statistical equilibrium with each other. This 
statistical equilibrium means that the turbulent motion in a given wave length 
receives its time -independent energy density from larger motions and discharges 
it to smaller scale motions. The general behavior can be observed qualitatively 
in many diverse turbulent flow phenomena. Consider as examples, the smoke 
from the end of a cigaret or the motion of the Gulf stream. In both cases the 
large scale eddies break up into smaller ones and these, in turn, disintegrate 
into still smaller motions. If energy continues to flow at a constant rate from 
the largest eddies, and to be dissipated at the same rate in the smallest eddies, 
the intermediate scale motions contain their fixed energy densities. Ultimately 
the energy reaches an extremely small scale, from which it deteriorates by 
shear viscous losses into the random molecular motion which we identify as 
heat. 

The band of turbulent velocity wave lengths through which the isotropic 
energy "cascades" from large scale to the smallest scale is called the 
"inertial sub-range". The smallest scale wave lengths within which viscous 
dissipation occurs is called the "dissipation" or "viscous sub-range " . 
Kolmogorov assumes that, within the inertial sub-range, the turbulence is 
statistically independent of the large scale sources from which it originally 
derived. In this inertial sub-range the turbulent power spectral density, 

(>d has been shown to be (Batchelor 1956) uniquely defined by the two 
parameters, the kinematic shear viscosity, p , and the rate of supply (or 
removal) of energy, c . The equation for the one dimensional PSD is 

$^(K) (2.1) 

where b = b (i^ , c) 

X = 2it/A= wave number of a turbulent 
velocity component 

A = effective wave length of a turbulent 
velocity component 

e = energy dissipation, per unit mass, per unit time 

l> = kinematic shear viscosity of water 
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The -5/3 rds power law has been amply verified for components of turbulent 
velocity in laboratory experiments and, to a lesser extent, at sea. 

We are interested in the turbulent velocity not only because it is a 
significant source of sound fluctuations in special situations such as turbu- 
lent channels but, more importantly, because it is a guide to the behavior of 
temperature variations . And temperature variations are generally the principal 
sources of sound fluctuations in the ocean. 

The temperature often acts as a convected, passive contaminent which 
closely follows the turbulent velocity behavior. The spectrum of temperature 
fluctuations has been studied for such a case (Batchelor 1959) and again the 
-5/3 law is predicted for the inertial subrange. The predictions have been 
largely verified (Stewart 1962, Grant 1968) by ocean measurements. 




Fig. 2 Temperature and velocity spectra at depth 27m in the open sea, from 

Grant (1968a). The predicted transition wave number, , is shown, 
as calculated from Eq . (2.8). 
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Fig. 2 plots the spectra (Grant 1968) of a component of velocity and tempera- 
ture measured at 27m depth, compared with the theories described in the 
references. The figure shows the extension of the Kolmogorov theory from 
the inertial subrange, where the -5/3 slope is confirmed, into the higher wave 
number dissipation sub-range. The boundary between the two subranges is often 
assumed to be at the Kolmogorov wave number. 



H = = (e 



( 2 . 2 ) 



Values of € , derived from measurements of the turbulence spectrum, are 
a f^nctioi^ of depth and have been found (Grant 1968) to be as high as 0.5 
cm sec st^an c^cean depth of 15 meters in well-mixed water and as small 
as 1.7 X 10 cm in quiet water below the mixed region at 213 meter depth. 

Looking again at Fig. 2 we see a rise in the temperature spectrum, also 
predicted ^y B^tchel^r (1959). At this depth, 2^m, ^ was determined to be 
5.2 X 10 cn^ /sec . Assuming v = 1.4 x 10 cm / sec , Eq . 2.2 yeilds 
X =6.6 cm . The rise in the temperature spectrum occurs between x and 

S -2 ° 

X "^2.4 X 10 X . At greater depths, the dissipation constant and the Kolmogorov 
wave number are generally smaller. 

Having identified the high wave number end of the isotropic -5/3 spectrum, 
we now seek to establish the low wave number limit. Unfortunately, at the low 
wave number end of the spectrum of temperature microstructure , the PSD is 
very dependent on the local conditions and history of the water mass, and only 
empirical results are available to guide us. 

A useful, but very crude, generalization that can provide some guidance 
is that the largest temperature wave length that could conceivably be 
isotropic at the depth H would be of extent 4H. The corresponding isotropic 
wave number would be 



m ■ 4H 



n 

2H 



(2.3) 



This would suggest that for Fig. 2 , where H = 2 7m, the lower limit of the 
-5/3 spectrum might be at 



H 

m 



5.8 X 10'^ cm'^ 



(or log 



10 



H =-3.2). 
m 



In f|ct, the data of Fig. 2 shojv that at values of x much smaller than 5.8 x 
10 (perhaps at x=“0.02 cm ) both the temperature and the velocity PSD 
break away from the -5/3 slope. This suggests that, realistically, the lower 
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end of the -5/3 spectrum to be expected at sea is somewhere between x and 

X . We will call that transition wave number, x . 
o t 

The appropriate variable to describe sound phase and amplitude fluctuations 
will turn out to be the speed of sound, and so it is the spectrum of fluctuations 
of the speed that we must seek. The connection is direct. For simplicity we 
will assume that the speed depends only on the temperature. We have previously 
seen that the fluctuations in temperature microstructure are much less than the 
average temperature. The fluctuations in the speed of sound are also a very 
small fraction of the average value. It therefore turns out to be most useful to 
define not only the index of refraction of the sound speed, at position R 

r\(R)= -- % ■ (2 . 4) 

<c(R)> 



but also the excess index of refraction, p (R) , which will be the particular 
parameter for all future discussions of fluctuations. 



^ (1). . W- < £ . (R)) . „ (R). i« I 

<C (R)) 



(2.5) 



The mean value, ( ) , is to be calculated over the time of the experiment. It 
is essential to realize that p (R) is sensitive to position K . However, for 
typographic convenience we will drop the notation for the explicit dependence 
on position so that henceforth we write p (R) =p . Further, we define 

/ \ j ( 2 . 6 ) 

(p) = o and o = \p ) 

P 

The point must also be made that if the time or duration of our averaging 
process is changed, the value of (p) and cr^ will change. We must immediately 
accept that these quantities must have their temporal limits restricted. A 
logical criterion for the restriction is that for the present we are concerned 
only with locally homogeneous regions of the ocean, only for times comparable 
to the duration of the acoustic experiment. We will be more precise about 
these limitations, shortly. 

Since the excess index, p, is proportional to the excess speed of sound, 
which in turn is assumed to be proportional to the excess temperature, AT, a 
spectrum of p can be expected to resemble the spectrum of AT. It will have an 
inertial subrange in which the log of the power spectral density, * (x) will have 

a slope of -5/3 when plotted against log x , as in Fig. 2. That isotropic inertial 
subrange will have an upper limit x* and a lower limit which we will call 
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The lower bound of the isotropic inertial subrange of temperature x , 
appears to be approximately a function of depth, H, alone if we define ^ 
it as 



K 



t 



0.5 



(h k ) 
o m 



0.5 (- 



v3 



)'/3 (- 



n 

2H 



(2.7) 



Fig. 3 is a graph of x^, calculated from Grant's (1968) values of c and H, as 
well as those of Nasmyth (1970). 




Fig. 3 Transition wave numbers, calculated from Eq. (2.7) using data of 
Grant (1968a) and Nasmyth (1970). The empirical fit to the data, 
given by Eq. (2.8) is shown. 



.The equation of the best simple 



approximation to this graph, is of the form 



where B = 9 . 0 m-1 



( 2 . 8 ) 



h = 40.0m 
H = depth , m . 

= lower wave number limit o£ isotropic 
temperature spectrum, m 
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For example, for the experiment at 27M depth, shown in Fig. 2, Eq. (2.8) 
allows us to diagnose that the medium was no longer isotropic for x<x = 

0.05 (logj^Q -1.3). Study of Fig. 2 verifies that especially the velocity, 
but also the temperature spectrum, breaks away from the -5/3 slope at 
log X —-1.3. We will use the criterion of Eq. (2.8) to guide us in our predic- 
tion of acoustic fluctuations due to microstructure. 



A graph of the spectrum of the expected excess index of refraction at sea 
is shown in Fig. 4. 




Fig. 4 Typical spectra ^ (x) and the approximate, simplified, version to be 
used here to calcMate spatial correlation. 



It is typified by: a large scale, (small wave number) non-isotropic, highly 
variable, section at x<x^ that depends on location and time of day: an 
isotropic universal subrange that follows the Kolmogorov Law 

-5/3 

4>^(x)=b(x) forx^^x^x*; 

-5 /3 

a universal dissipation subrange with $ (x) >bx for x* ^x^x , and a 

rapid exponential drop-off for x>x . will generalize and use thus spectrum 
in the next section. The data avaiPable show that x «x. «x^. 

o L m 

2 . 2 The Spatial Correlation and Structure Function 



The alternative to the spectral description of fluctuations of temperature 
or index of refraction is either the correlation function or the structure function. 
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which is simply derived from the spatial correlation. These functions turn out 
to be the most direct way to calculate the sound fluctuations , as we will see 
in Section 3 . 

One of the earliest determinations (Liebermann 1951) of the spatial 
correlation of temperature fluctuations was conducted by mounting a therm- 
istor on a submarine operating at a depth of approximately 50m. The 
standard deviation of the temperature fluctuations was found to be approxi- 
mately 0.4C at this depth (time of day, weather conditions unreported). 

The spatial correlation could be fitted approximately by the exponential law 
(shown by solid line in Fig. 5) 

(5) s ^ 

where a 



= e 



-§/a 



= 60 . cm 



(2.9) 



which thereby defined a correlation length, sometimes called a "patch size", 
of 60 cm . 




Fig. 5 Temperature correlation function from Liebermann (1951) with data 
(circles), empirical fits to exponential (solid line), and Gaussian 
(dashed line, correlption functions. The correlation based on Fig. 4 
with = 0.24 cm is shown by dotted line, identified as "Kolmogorov" 



Others have used the Gaussian correlation function 



Ct (4) 



e-U/a)' 



( 2 . 10 ) 



to describe microstructure fluctuation. A Gaussian correlation function 
(dashed line) has been drawn in Fig. 5 for comparison with the exponential 
functions and the data (circles). 
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Bo^h functions define the same correlation length, since they both reach 
e at ^ = a. The correlation function of the index of refraction will equal 
that of the temperature, if other parameters affecting the speed of sound 
have negligible effect, and if the fluctuations are small. Then 

Cj («) (?) 

It is not always possible to determine the correlation function C,j, U) or 
C (4) by measuring it directly. In any event it would be a boon if we could 
u^ the known universal characteristics of ocean microstructure to deduce 
the spatial correlation function in the real world outside of the laboratory. 
Although this turns out to be not completely feasible, there are certain 
generalizations about the spatial correlation function that can be deduced; one 
of these is its form at small values of 

The spatial correlation function of the excess index of refraction can be 
obtained from the wave number spectrum by the Fourier transform theorem. In 
three dimensions this relation is 

= —1— Iff dH (2-11) 

a ^ 

where T is the vector spatial lag. 

If we now assume that the changes of the index of refraction are independent of 
the direction being studied we can simplify the integration. In spherical 
coordinates 



—i# -4 2 . 

ht. = and dn = 2nH sinO- dn d 1 ^ 

Then C (r) - j" k2 

a 

M- 

and the integration yields 

^ (H) JX 

We now define the one dimensional spectral function, $ (x) for the isotropic 

medium ^ 

(h) dH = 4tt (h) dH = $^(h) dH and change (2.12) 

notation, r ^ , for consistency with previous work, so that 
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00 

(5) - Jo (K) dH. (2 . 13) 

^ may be in any assigned direction for the isotropic medium. 



In order to integrate (2.13) we start with the spectral conclusions of the 
previous section (Fig. 4). We are prepared to make four simplifying assump- 
tions in order to proceed: 



1 ) $ (k) =$ +s(h-h) for H < H 



2 ) 

3) 

4) 



^mr^o 

where s is the slope constant= *' 

= peak isotropic PSD; ^ (o) 



$ (w) 

M- 



constant = 



$ (h) = $ 

|X ra ' 

$ (h) = o 

It 



„ for ^ ^4. 

m m t 

for < K < H 
t o 



for 



H > H 



(2.14) 

Source Subrange 

Transition Subrange 
Inertial Subrange 
Dissipation Subrange 



The approximate spectrum that follows these assumptions is shown by the 
dashed line of Fig. 4. 

The effect of the simplifying assumptions is to replace the highly variable 
anisotropic range and the fixed, but complicated, isotropic range by a fixed, 
easily-integrable , spectrum. This will permit us to assess the relative impor- 
tance of each of the sub-ranges and to determine what universal generalization 

can be made, (if any) about C (^) at sea. 

b 

Using the assumed spectrum of Eq. (2.14) it ^s possible to immediately 
determine the excess local index of refraction, a in terms of the peak spectral 
density, $ , and the constants and x^. Recahl'ling that the mean squared value 
of any fluctuating quantity may b^ obtained by integrating the spectral density over 
all values of the appropriate parameter (frequency or wave number) we find 









§ (h) dH 

O p. 



= $ 



ra 






1 + 



(X - 



H )] dH + J dH + J, (— ^) dHj 



m 



m 



K 



m 



H. 



H 



Since x x, and $ < $ , the result of the integration simplifies to 

m t o m 



2 

^11 = 5/2 K $ 
t m 



(2.15) 
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It is interesting to determine the parts of the spectrum most responsible 
for this variance of the fluctuations . The terms in the integration were 



1/2 X $ source subrange 

m m 

X, $ transition subrange 

t m 

3/2 x^ $ inertial subrange 

t m 



Since x « ^ 4 . conclude that the contributions to a come almost equally 

from the transition and the inertial subranges. Furthen^ore, our simplifying 

assumption of the particular form of the source subrange can have little effect 

on o because that contribution is proportional to x , and x «x . 

|j. m m t 

A second general conclusion can be reached if we restrict our interest to 
very small values of the spatial lag, 4 • Then we can substitute the expansion 
for the sin and integrate Eq. (2.13) with ease 



C 

M- 



(§) 



^2 






Jo 



$ (h) dx - 






il 



.H 

m 









L 

$ 

m. m 



for K § « 1 
o 

(2.13a) 

[§ -f m (k - H )] dK 
dK 



+ 




H 



dH 



The first integration yields (fj, ), as we saw above, so that the first term 
is unity. The integration over the three terms in the curly bracket is also simple 
so that 



C 

M- 



(?) 



1 



-2 , 5/3 A/3 

g Ht Ho 




(2.16) 



and with the aid of (2.15) 



C (§) = 1 - 0.05 k for H I «1. (2.17) 

to o 



The significance of this result appears _if we consider any specific depth ^ say^^ 
50m. Then (2.3) gives x = tt/IO cm , (2.8) predicts x = 2.4 x 10 cm 
and therefore, at 50m dep?h, we can expect 

(?) = 1 - O.I 168 §2 
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For small displacements this is an expansion of the Gaussian form 



C 

IJ* 




(5/a)^; 5 « a 



where ^ ^ 



We therefore conclude that at a depth of 50m the Gaussian form is 
appropriate only for very small spatial lags, ^<sc0.68 cm. The conclusion is 
not significantly different at other depths. 

We are now prepared to develop the universal correlation function based 
on the spectrum of (2.15). The correlation function will be "universal" only 
in so far as it is calculable from the universal part of the spectrum. To be 
safe, we will calculate the correlation only for ir /x^. 

For this restricted part of the correlation function it turns out that the universal 
Kolmogorov subrange is dominant so that very useful results can be obtained 
with a minimum of mathematical calisthenics. 



Starting again with Eq. (2.13) we split the integration into the contributions 
for the source, transition, and inertial subranges. We have assumed zero spec- 
tral density in the dissipation subrange so that integral is zero. 



a C (§) 
ii IJ- 



' C"- „ + s (H - K ):a« + J " + J 5/3 

K = o TO ' m^-* K m X§ to'h 

m 



TO^H 



The coordinate for integration is now changed to |3 = x ^ 



Hi 



M- IJ- " " § 3 = o'- P TO I m 5 r=k E P P=H I m 



■sin P- 






The approximations are now made for J d/9= j9- r^rr. +■ 



ll 



in the source 



3.3! 5*5 1 

integral and in the transition integral. Because of the rapid convergence due 
to the j8 8/3 in the denominator, the third integral is comfortably approximated 
by using the sine expansion from zero to ^72 and again from ir/2 to it. Using 
(2.15) the result of these integrations, valid for ^ with about 5% error, is 



2/3 



C (^) = 1.0 - 0.48(x^^) ^ +0.28(x^4) 



+ 6 . 7 X 



10 



(2.18) 
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The curve is shown in Fig. 6 
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Fig. 6 Spatial correlation function C (x derived from spectrum of Fig 4. 

M" t 



It turns out that the source subrange makes essentially no contribution to 
the correlation function for 4 -it/x^; its effect would be for greater spatial lags. 
The transition and inertial subranges, together, are totally responsible for 
C ( 4 ) for ^ ^tt/x . The inertial subrange is particularly dominant for small 
si^tial lags (the 2/3 power term). 

Assuming that is given by Eq. (2.8), it is possible to calculate C ( 4 ) 
for any depth. The calculation has been made for the depth ^f LiebermanH's 
experiment (50m), for which Eq. (2.8) gives x^ = 0.024 cm , and the predicted 
correlation is plotted as the dotted line in Fig. 5. The agreement with the 
experimental data is remarkably good, up to the limiting spatial lag, 

Recapitulating, the universal spatial correlation function that has been 
derived from the assumed isotropic ocean spectrum has a Gaussian correlation 
function from ^ = 0 to approximately 4^2'nr/x . It then drops off with a slope 
of close to -2/3 under the control of the inePtial subrange. This slope decreases 
to close to zero at the limit of our study, ir/)^. Experimental studies show 
that, beyond the correlation value will continue to decrease slowly, 

begin to oscillate about this drop-off with a spatial periodicity and a specific 
character that depends on the precise form of the spectrum in the transition and 
the source subranges. 
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3. SOUND PHASE AND AMPLITUDE FLUCTUATIONS 



3 . 1 Phase Fluctuations 

Let us begin our study of the variable propagation of sound in the ocean 
by considering plane wave sound phase fluctuations along a ray path in a 
statistically homogeneous, isotropic, medium. We will assume that the 
sound frequency is high enough so that the wave length is much smaller that 
the correlation length of the microstructure, that is ka »1. Again a is th^ 
correlation length, the distance in which the spatial correlation drops to e 
We will first consider the ray approximation, and relatively small propagation 
path, so that ka »R/a, where R= path length. 

The dual condition is best written in terms of the wave parameter. 




For this ray development (Bergmann 1946) we require D «1. The time to 
travel the distance, R, in the x direction is given by 



where n = c/c = local index of refraction of the medium, T = clock time and 
c^ = (c) . BotS c and n are functions of position and time of the experiment. 

The particular duration required for this traverse will depend on the type 
and arrangement of inhomogeneities at the instant of the sound travel. The 
time of travel is assumed to be short compared to the time for the medium to 
rearrange itself. That is, the medium is assumed to be "frozen" during the 
sound traverse , and to change into another configuration during the time 
between traverses. 

Our interest is in the duration of the travel, averaged over either many 
statistically equivalent inhomogeneous paths, or averaged over many traverses 
along a time-varying inhomogeneous path. We assume that these two averaging 
processes are equivalent, that is, we assume that the medium is ergodic. 




dx 




o 



o 



because the average 



index, <n> , is unity. 
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It is the deviation from the mean that we must calculate: 



At = t - (t) :r J n (x,T) dx - J dx 

CO CO 

o o 

1 

= — p.(x,T)dx 
where (x,T) = n (x,T) -1 



From this we get the square of the deviation from the mean (At)^ 
^ tidx which can be written (At)^ ~~c^ Jo i*o 



dXfdXz 



2 1 R R 

The mean square deviation from the mean is <At > = — - f dx^ f Ai(x, ,T) 

rT,\\, c2Jol*JOl 



^ (x^,!)) dx^ 



The integrand can be written in terms of the spatial correlation of the excess 



index of refraction C ( 4 ): 

- M' 



<p,(Xj^,T)p,(x2,T)> = (p,^) C (x2~Xj^) = 



<At^> = C 



In these terms 



c 2 
o 



O [X 



(X2-Xi>dx2 



(3.2) 

(3.3) 



We seek the mean squared phase fluctuation which is calculated by 
substituting <p= wAt where (pis the instantaneous deviation of the phase from 
the value due to the mean speed, c . The units of <p are radians. Since the 
medium is presumed to be statistic^ly homogeneous, it is the separation, 

^ = 1^2 “ ^1 I' significant not the particular end points, x^ or x^. 

Therefore, The mean squared phase fluctuation is 



<(p^> =a^k^.f^ C ( 4 ) d4 D«1 



(3.4) 



2 2 2 2 
Since <cp > = (cp) + Var <pand <(p) = o, by definition, <(p )= Var cp= . (3.5) 



We have only to know the spatial correlation of the index of refraction, from 
the source to the receiver, in order to carry out the integration, and obtain 
the variance of <p . 

2 2 

It should be pointed out that the quantity a = <|i > must be known at the 
time, and location of experimei^t. The signific^ce of this comment is that, in 
the ocean, the evaluation of a can be strongly dependent on position and the 
range at a given time, or a function of the time and duration of measurement. 
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For example during passage of an internal wave, at positions in the thermo- 
cline there will be large variations of the temperature and the speed of sound, 
over times of the order of several minutes. 

It is also convenient, sometimes, to assume that the spatial correlation 
function drops to zero after a few correlation lengths in either direction from 
any reference point. Then for the sound path of length R »a we can extend 
the integration to infinity in both directions , or double the integration from 
zero to infinity, without affecting the total integration, 

00 

a^= <tp^> = 2 k^Rf C H) dP D« 1 (3.6) 

(P U V 



Given a correlation function that drops rapidly and monotonically to zero, 

Eq. (3.6) is more readily evaluated for the mean squared phase fluctuation 
for the high frequency condition than (3.4). 

When the path is long, so that phase interferences occur, a wave analysis 
is necessary. We will not go through the derivation for the phase fluctuations 
through the same type of medium at this other extreme of propagation parameters , 
but simply quote the result (Chernov 1967), applicable for low frequencies or 
very large path lengths. 

k^ Rf”c (4) dF D»1 (3.7) 

<p p Jo p ’ ^ 



The point to be observed at this time is that, regardless of the value of 
the wave parameter, D, the variance of the fluctuations of phase are proportional 
to the length of the sound path, the sound frequency gquared, and the variance 
of the fluctuation of the excess index of refraction, a / along the path. Further- 
more, regardless of the form of the spatial corre latioi^ function , there is only 
the factor, 2, relating the mean square phase fluctuation obtained for very 
large or very small values of the wave parameter, D. This factor will not be 
altered if we change the form of the spatial correlation function, either for 
convenience in the evaluation of the integrals , or as we seek to describe 
the medium more accurately. 



Another point must be made at this time. Mintzer (1954) has shown that 
solutions such as we are considering can be valid only for single scattering 
in a weakly in homogeneous medium. That is, although we have ka »1 as a 
requirement for the relative size of the "patches", and kR as a requirement 
for far field propagation, the fluctuations of the medium, o must be so small that 



k 



2 



2 „ 



« 1 



In practice this c^onditic^ is not difficult to fulfill; field measurements 
generally show <10 
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The Gaussian correlation function 



C 



H) =e 



-(4/a)^ 



(3.8) 



is attractive because of the simplicity in evaluating the integrals of Eqs . (3.6) 
and (3.7). For the Gaussian correlation function the phase fluctuations are: 



a^ = ((p^) = yiTa^k^ Ra D«1 

<p b 



(3.9) 



and 



<P 



/ 2. yiT 

= <<p > =— 



2 ,,2 



R a 



D» 1 



(3.10) 



It has been shown (Krasilnikov 1956, and Chernov 1967) that, for a 
Gaussian correlation function, the variance of phase fluctuations over all 
values of D is described by 



R a ( 1 + ^ tan’^D) (3.11) 

0-2 u 



3 . 2 Amplitude Fluctuations 



The amplitude fluctuations are often put in terms of the ratio of the 
change of the pressure amplitude to the original amplitude, 

V = ( ) (3.12) 

P 

o 



where P is the sound amplitude at x = R and P is the amplitude of the plane 
wave at the sLarting position x = o. In terms of this parameter it turns out 
(Bergmann 1946, Mintzer II 1953, Chernov 1967) that the mean squared 
fluctuation of amplitude in the high frequency extreme is 

(3.13) 

and in the wave interference region (Mintzer I 1953) it is simply 

^ <x^ = a^ pf C (^) d^ D»1 (3.14) 

A [1 Jo p ^ ^ 

where o is the fractional standard deviation of the pressure amplitude 

/A 




2 , 2 . 

cr=(x > = 



_kL 



R 






2 

for D - 4RAa « 1 



a 



A 



<( 



P-Pq 

Po 




(3.14a) 
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Comparison of expressions for amplitude fluctuations with those for phase 
fluctuations is revealing. For example, for the high frequency situation D «1 
the mean squared amplitude fluctuations depend on an integration of the fourth 
derivative of the correlation function in the transverse directions, y and z, 
rather than a straightforward integration over the correlation function in the x 
direction as in Eq . (3.8) for the phase shift. This means that the amplitude 
fluctuation is caused by the curvature (really double curvature) of the correla- 
tion of the refractive index, transverse to the sound path. The dimpled medium 
in the case of amplitude fluctuations therefore behaves as a series of lenses 
coverging or diverging acoustic energy along the sound path (notice the evalua- 
tion at zero displacement in the y and z direction, that is , along the beam). 

At the other extreme of propagation parameters, however, when the range 
is larger, or the frequency is lower, such that R/a » ka (D »1), the sound 
field fluctuations are caused by field distortions due to a large number of 
inhomogeneities along the path in the medium; these distortions interfere, 
constructively and destructively along the entire path and this interference 
effect is common to both phase and amplitude. For this condition the mean 
squared fluctuations of phase and amplitude are identical. Both are linearly 
dependent on: a) the distance of propagation, b) the frequency squared, and 
c) the integrated spatial correlation function. 

It must be pointed out that all of the derivations to this point have been 
based on the assumption of single scattering by weak fluctuations. The 
consequent solutions are valid only in so far that very long paths and very 
large fluctuations are not allowed. In practical propagation over very long 
paths the magnitude of the phase and amplitude fluctuations cannot increase 
with R without limit, or they would reach the absurd conditon of fluctuations 
that are greater than the quantity itself. In fact, "saturation" must take 
place and the increase of the fluctuation with range must reach a limit that 
will be indicated, shortly. 

If the Gaussian (^orrelation function of the index is assumed (Krasilnikov 
195 6) the variance, a , is: 



<X^> = k^ R a (l-^tan“^ D)] 

<X^> = (-^) (R/a)^ D « 1 

% k^ D » I 



(3.15) 

(3.16) 

(3.17) 
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3 . 3 Small Amplitude Fluctuations: Laboratory Experiments 



The dependence of small values of on the parameters of the medium and 
the frequency has been put on a firm experimental basis in a beautiful set of 
laboratory studies (Stone and Mintzer 1962) in which the microstructure of the 
index of refraction was obtained by heating the tank of water through which 
the sound beam was propagated. First the rms value and the spatial correlation 
function of the temperature ^icrostructure were obtained by direct measurement. 
From these data the value o , and the correlation function, Cp ( 4 )/ of the 
excess index of refraction \lvere calculated. The spatial correlation function was 
fitted to the Gaussian form, and the correlation length, a, was calculated from 
Eq. (3.8). The measured temperature microstructure, fitted to the Gaussian 
correlation function showed a correlation length, a = 3.5 cm. 

Additional comparison showed that the rms excess index of refraction, ^ 
deduced acoustically, was approximately 25% less than the value 1.6 x 10 , 

calculated directly from temperature measurements. One flaw was revealed 
in the experiment; the acoustically determined value of the assumed Gaussian 
correlation length was only about one-third of the value (3.5 cm) directly 
measured by thermistors. This discrepancy was possibly due to the error of 
assuming that the correlation function was Gaussian. As shown in section 2, 
the Gaussian form is appropriate only for very small displacements. A 
combination of Gaussian and Kolmogorov correlation functions would probably 
have been a better description of the medium for insertion into Eqs . (3.13) or 
(3.14). 

3 . 4 Large Amplitude Fluctuations - Large Ranges 



So far we have considered small fluctuations for two cases: 

1) Change of sound pressure amplitude, P, in the presence of zero 
(or negligible) phase shift. In section 3.2 we gave expressions for 



2 _ P-P 2 

""a <ej-^) > • 



The type of pressure field at x = R that results in is 



o 

described by the phase diagram for”P"in Fig7a where P represents two 

possible scattered pressures at time, t position x. It is the vector addition 

of various values of P to the undisturbed field, A (x), that results in the 

statistically changing resultant pressure am^itude, P, at angle <p , with the 

original wave. This P is used to calculate . It is to be noted that the 

random amplitude P has components, a , in phase with the original wave, 

A, and b , 90 out-of-phase with A. It^is a rather than b that is effective 
s s s 

in changing the value of P. 



2) Change of sound phase due to varying index of refraction. In the 

derivations of section (3.1) leading to expressions for q it was assumed that 

amplitude fluctuations have a negligible effect on phase fluctuations. From 

Fig. 7 a it is seen that the out-of-phase component of scatter b , rather than 

s 

, is the principal source of phase changes. 
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(a) SMALL FLUCTUATIONS 




(b) LARGE FLUCTUATIONS 

Fig. 7 

Two realizations of scattered pressure, P added to undisturbed pressure A. 

P is composed of the in-phase part, a and the out-of-phase part, b . For 
small fluctuations Fig. 7a the pres sure ^amplitude P, is principally a mnction 
of A and a only, and the phase (pis principally determined by A and b^. For 
large fluctuations Fig. 7b both P and cp depend on A, a and b . 



We now consider large fluctuations for which the situation of Fig. 7b is 
more appropriate. We will no longer consider that the phase fluctuations have 
negligible effect on the mean square amplitude or that the amplitude fluctuations 
have negligible effect on the mean square phase. 

The complete picture at x = R, which we have represented in Fig. 7 can be 
described analytically by 

p(x,t) = A(x) cos (cot - kx) + a (x,t) cos (cot - kx) + b (x,t) sin (05t - kx) 

(3.18) 
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where A (x) is the magnitude of the pressure of the remaining unscattered field 
at x; 

a (x,t) = time varying magnitude of the scattered pressure which is in phase 
with the original wave at x. We assume that values of a are Gaussian distributed 
and that <a > = 0; ^ 

b (x,t) = time varying magnitude of the scattered pressure which is 90 out-of- 

p^ase with the original wave at x. We assume b is Gaussian distributed and 

<b > = 0. b ^A. ^ 

s s 

Equation (3.18) resolves the fluctuation scattered pressure, P of Fig. 7 into 
components in-phase with the depleted original wave A (x) ancf out-of-phase 
with it. A typical situation in Fig. 7 shows that 



t .2 2 ^2 

P = a + b 
s s s 



To get the typical amplitude, P, of the total wave at x as a function of A, 

(3.19) 



a , b , Eq. (3.18) is rearranged to 
s s 



P = (A + a ) cos (tot - kx) + b sin (wt - kx) 
s s 



= P cos (u)t - kx - «) 

where P = pressure amplitude at x 

r/. ^ 2^1/2 

= [(A + a ) +b ] ^ 
s s 

$= phase angle with respect to original 

-1 .^s . 

wave at X = tan (. . "•") 

'A+a 

s 



It is useful to expand ^ and $ for small fluctuations a^ 
Using the binomial theorem for P 



(3.20) 



(3.21) 



b « A. 
s 



P = A(1 + 
=“ A(1 + 



2a< 



. £al+bsi_) = 
+ 2F ’ 



A+a 



( 3 . 22 ) 



and the small angle approximation for $ 



* -1 / bs X 



(3.23) 
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It is therefore clear that for small fluctuations the o’, of Eq. (3.12) and 
(3.17) can be interpreted from the point of vie'w of (3.18]Tand (3.22) 



^ = -P-Pp- 
Po 



A + ag - A a 



... 2 ...» Var a <5 

so that = Var(x) = 2 -^- 

o 



(3.24) 



As we anticipated from our study of Fig. 7a, the value of cr^ is therefore 
relatively independent of the out-of-phase components of scatter, b , because 
their effect shows up to be relatively small in the small fluctuation approximation 
of (3.22). 



Likewise, for phase fluctuations the o of Eq. (3.4) - (3. 11) can be 
understood from (3.18) and (3.23) as ^ 



a = 



Var be 



A 



r 



(3.25) 



which shows that the value of o is virtually independent of the in-phase 
variations, a^ , in the small fluctuation approximation of (3.23). 

For large fluctuations we follow the method of Brownlee (1973). First 
we recognize that the diversion of acoustical energy into scattered energy will 
cause the average in-phase pressure amplitude to reduce from the value A 
at X = o to A at x; we will assume that the rate of attenuation from the 
original magnitude is proportional to the distance traveled. 



dA 

A 



- 0 / dx 

s 

where a = attenuation of the original 

^ wave due to scattered energy. 



Oi R 

Integrating for a range R, A = A e s where A = P = pressure (3.26) 
amplitude of (unscattered) wave at x=o. 

The average intensity of the part of the wave in step with the original wave 
will therefore be proportional to 



A 



2 




-2a R 
e s 



(3.27) 



We will assume conservation of energy for the total wave system; that is , we 
assume that whatever energy is taken from the original wave goes into the 
scattered wave and that this is principally forward scattered (Skudrzyk 1961). 
This means that the intensity is conserved, when averaged over any plane 
perpendicular to the direction of propagation . 
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From (3.20) 



<pQ^> = ( P^> 

<Aj^ > = <(A+a )^ + b 

= <A^> + <2a^ A> + <a ^> + (b 

"o = ^ <-s'> ^ <^'> ' 



To determine the attenuation rate, a we use (3.27) and (3.28) and form 



An^-A^ ^ (as^> + (bs^ > _ , .-20! R 

2 2 
A 

o o 



(3.29) 



The evaluation can be made in the small fluctuation range where the 
ingredients of (3.29) are well known. There, using the expansion of the 
exponential for small argument and (3.24) and (3.2 5) 



= 20! R 
A cp s 



(3.30) 



The condition for long range is of particular interest. 
Then we have, from (3.7) and (3.14) 



a 



2 

A 




» 1 



(3.31) 



so that the attenuation constant of the original beam due to energy going into 
fluctuations , is 






s 




^ D » 1 



(3.32) 



In particular, when the scatter is due to a medium with a Gaussian correlation 
function, C ( 4 )# we have, using (3.17), 



OL 

s 



2 J- 



(3.33) 
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We can now proceed to consider large ranges and large fluctuations. 

The quantity observed in any experiment is the pressure amplitude P. Eq. (3.20) 
shows P to be a function of all three quantities: The steady coherent component, 

A, and the fluctuating random components, a and b , The conventional way to 
express the amplitude variations is in the ratio ^ 



(CAV) = 



<P^> ->)' 

~2 



<P>‘ 



(3.34) 



where CAV= coefficient of amplitude fluctuations. It might be noted that for 
small fluctuations we have simply CAV “ c^» Eq. (3.14a). Here, for large 
fluctuations we find, from the conservation of energy assumption that 

<P^>= <P^^> = A^^ (3.35) 

2 

But the calculation of (P> involves a more involved evaluation of the mean 
value of P given by Eq. (3.20) when both a and b vary randomly, each with an 
assumed Gaussian probability density function. ^ 

The mean value of the sum has been found by Brownlee (1973) and is 




where S = 



-e ^®sR]^[exp[-2/ (e^^^R-1)]} ^ 



1 + 3b + 



(3.5)b2 

(21)" 



+ 



(3.5.7) b3 

(Jl7~ 



z 



b . -l)'^ 



(3.36) 



2 2 

When this is now squared and combined with P in Eq. (3.34) the A cancels and 
the result is expressible solely in terms of the attenuation constant, a / ^nd the 
length of path, R. 



(CAV)^ = -1 +-1 exp ( dV V ^ (-2a R)][i:^]r^ (3.37) 

7 T exp(ZQ! R)-l S 

This messy expression has a simple maximum value for large fluctuation and 
large ranges; we get 

(CAV)2 . -I + - 0.273 £cr2a^R»l 

or (CAV) = 0.522 
max 

The complete behavior of the coefficient of amplitude variation is shown in 
Fig. 8. 
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Fig. 8 Coefficient of amplitude variation as a function of the standard deviation of amplitude fluctuations 
for small fluctuations, o^. The data of Shvachko (1967) at 25. kHz are shown. The fluctuations 
are assumed to be due to microstructure only. 



which is the graph of the square root of (3.37). We observe that the fluctuations 
(CAV) = 0^ for small values of that is, for small fluctuations. The linear 

dependence of (CAV) on range continues until approximately = 0.35 after which 
growth slows. (CAV) reaches 90% of its ultimate value at approximately 0 = 

0 . 6 . ^ 

2 

Knowing the microstructure constants, a and a. Fig. 8 and Eq. (3.37) 
may be used to determine the range at which 4he CAV due to microstructure can 
be expected to reach any value up to the maximum CAV = 0.522. 



3 . 5 Sound Fluctuations at Sea 



In section 2 we found that the beautiful, and easily manipulated, Gaussian 
correlation function is precisely correct only for very small spatial lags at sea. 

In the previous section we saw that the Gaussian correlation function may yield 
the wrong answers for the coefficient of amplitude variation, even when very 
carefully tested in the laboratory. It is time to turn to the expression for the 
limited microstructure correlation function developed in section 2 for an assumed, 
reasonable, ocean spectrum , Eq. (2.18). 



The standard deviation of the phase fluctuations and the coefficient of 
amplitude variation for all values of the wave parameter, depend on the inte- 
gration of the spatial correlation function over all values of ^from ^ = 0 to 
infinity. Let us look realistically at the upper limit of the integrations such 
as in Eq. (3.7) and (3.14) which yield a and 0^. However, although the 
correlation, defined by ^ 



C^(S) 



<u(x) p(x+(i)) 

a 2 
d 



can be determined for any spatial lag, the value of C (O quite 

irrelevent to a calculation of the sound fluctuations f^r a path of extend, R. 
Therefore practical acoustical application of integrals of C (4)^4 
satisfied by an upper integration limit, R, rather than <» . 

The microstructure study of section 2.2 leads us to have confidence in the 
form of C (4) only for spatial lags over which the Kolmogorov spectrum is 
dominant.^ We have set these limits as 0 In order to evaluate the 

remainder of the integral consider Fig. 9 which is based on temperature meas- 
urements by Seymour (1971) and transformation to spatial coordinates by 
Haley (1972). The experiment was performed at a depth of 7.3 m, using the 
stable NUC Oceanographic Tower, 1 mile off San Diego. At this depth, 

Fig. 3 indicates a micro structure with transition wave number, x ^ 0.075 cm 

Equation (2.18) then gives the Kolomogrov-dominated spatial correlation function 
plotted as the dotted line curve to its expected limit ir/x^ = 42 cm. The 
agreement with the experimental data is excellent until ^==25-30 cm, after which 
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Fig. 9 Spatial correlation function of temperature microstructure at a depth of 7.3m. The solid line is 
the experimental measurement. The dotted line is the truncated universal correlation function 
based on the Kolmogorov spectrum from Fig. 6 with ^ = 0.075 cm . The dashed line is the 
Gaussian correlation function with the same correlation length, 25. cm, as the Kolmogorov. 



it levels off while the actual correlation function shows a continuing decrease 
and then an oscillation with small amplitude around the zero correlation axis. 



In order to obtain the integral of C ( 4 ) for a propagation experiment in such 
a sea we have two choices: One solution would be to determine C^(4) and 
C ( 4 ) out to the range of the acoustic experiment; this would not only be a 
formidable task but unrealistic for any long range experiment. The second, more 
expedient, choice is to accept a judicious combination of theory and empiricism. 
We propose to use the theoretically based microstructure of Fig. 6 to determine 
the correlation length, a, and then to use the Gaussian correlation function to 
approximate the true correlation curve. Judging by Fig. 9, this will result in 
slightly too large a contribution for X C (^)d^ between 4 ~ 4“ 0.5M. We 

will assume that subsequent positive aHd negative parts of C (4) will cancel 
each other and contribute nothing to Xc ( 4)^4 from 4 0 • 5 to 4 = R. In the 
example being discussed, the value of ttie correlation length, a = 25 cm. 



Having found a suitable way to determine the equivalent correlation length, 
all of our algebra for the sound fluctuations due to a Gaussian correlation func- 
tion is assumed to be applicable, even though the true correlation function is 
by no means Gaussian. 

The general method for the determination of the correlatiop length is found 
by referring to Fig. 6. The correlation function is down to e at 



x^4=x^a = 1.8 



(3.39) 



therefore for our approximation 



1.8 




Furthermore, using the empirical approximation to the transition wave number 
that we developed from Fig. 3, Eq. (2.8), we have the generalization for the 
upper ocean. 



a = 0.2 e 



H/40 



where H = depth, m (3.40) 

a = temperature (or index) correlation distance in m. 



It is evident from the data scatter of Fig. 3, from which the constant 
B and H are determined, that (3.40) can be only a guideline for the prediction 
of the correlation distance. 
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Although a often decreases with increasing depth, a general rule describing 
(j as a functidn of depth is not possible at this time. In fact, when such a rela- 
tion is found there will certainly be other important parameters in the equation, 
such as time of day, to which cr^ is particularly sensitive at shallow depths. 

When a is measured at the place and time of an acoustics experiment, 

Eq. (3.40) M'ovides a value of correlation length, a, that completes the needs 
for evaluating ql or a . For example, Shvachko (1967) has made measurements 
of a and CAV ardepMs 20, 35, 40M using sound of frequency 25 kHz. If we 
tak^Schvachko’s value of cr , R/ k, and use his depth in Eq. (3.40) to obtain 
a, we can calculate cr^ from*^(3.17) and predict his CAV from (3.37) or the graph 
Fig. 8. This has been done in the following way: 

ILLUSTRATIVE EXAMPLE 



Shvachko (1967) has performed an acoustic fluctuation experiment at 
frequency 25 kHz, range 20 m, in the sea of Norway. The medium ^ad a 
standard deviation of the excess index of refraction q = 17.8 x 10 . We are 



to calculate the predicted q . 



1-1 



-1 



A^(^/^0pth 20 m reading from Fig. 3 = 0.055 cm or using (2 . 8) x ^ = 

9 . Oe = 5. 45 m 

1 8 

therefore from (3.39) a = — ■ — = 0.33 m 

t 

or directly from (3.40) a = 0.2 = 0.33 m 

^ /l^ 7 (7.S Y 



0) 2 (25.x 10 ) ... - -1 

= - — = = 104.6m 

1.5 X 10'^ 



calculate D= ~ 

ka^ = (104.6)(.33^) 

therefore, using the Gaussian assumption (3.17) for D »1 

^ o ^k^Ra 
A 2 u 

= (0.89) (17.8 X 10'^)^(104.6)*(240.)(0.33) 
= 0.0244 
= 0.16 



This is to be compared with Shvachko's CAV = 0.17 
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Since the a, is small, we are in the region of the curve (Fig. 8) where we 
should expect CAV = o^. Therefore the small discrepancy between our predicted 
and the experimenter's measured value of CAV is a good confirmation of our 
generalized theory based on the truncated universal correlation function of the 
temperature microstructure. All of the other data points except one, from the 
Shvachko experiment, provide similar confirmation as shown in Table 3.1 and 
Fig. 8. 



Location, Time 


Depth m . 


R.m 


q X 10^^ 


CAV 


Theory 
q^ (-CAV) 


N.W. Atlantic, May 1962 


40 


200 


13.7 


0.095 


. 136 


ft 


40 


665 


25.4 


0.264 


.46 


II 


40 


1150 


4.7 


0.141 


. 113 


N.E. Atlantic June 1962 


35 


940 


2.7 


0.071 


.055 


II 


35 


600 


5.4 


0.100 


.087 


It 


35 


480 


11.0 


0.158 


0.16 


Sea of Norway, July 1962 


20 


160 


17.8 


0.173 


0.16 


II 


20 


160 


25.2 


0.173 


0.18 



Table 3 . 1 Comparison of experimental data of Shvachko (1967) with predictions 
based on the microstructure theory Eq. (3.40 and 3.17) of section 3. Sound 
frequency, 25 kHz. The points are plotted in Fig. 8 where is the theoretical 
value and CAV is the experimental value. The symbols in the first column are 
those used on the figure. 

3 . 6 The Limits of Predictions 



The calculation of amplitude fluctuations by the techniques of this report 
has been constrained by several approximations and assumptions. The ex-post- 
facto type of comparison with the experiments of Liebermann and Seymour- 
Haley show that our truncated universal temperature correlation function based 
on the universality of local microstructure in the upper ocean is not far from 
the mark, and is usable as a function of depth, alone. 

The highly variable values of o in the ocean preclude our use of the 
equations for predicting sound amplitude or phase variations unless we have 
in-situ data of q at the time of the experiment. However, if q is available 
as in the Shvachhco experiments , our predictions of and CAV^^are quite 
acceptable . 

We have continually emphasized that our predictions are limited to the 
upper ocean. By this we mean approximately the upper 100 m in temperate or 
tropical waters. The reason for this restriction is that the recent data of 
Nasmyth (1972) for depths below 100 m do not fit our Fig. 3 which we use 
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to determine and correlation distance, a. The reason is clear, at greater 
depths in the temperate oceans, and even at shallow depths in the more stable 
polar regions , the common existence of layers and thin sheets of highly 
distinctive water masses contradicts our requirement for locally homogeneous , 
isotropic water. Under these conditions it is not yet possible to generalize in 
the sense of this section. 

A further limitation to our development is the continuing restriction, 

ka » 1, on which the derivations for 0 ^ and a were based. At shallow depths, 

say H=10 m, the correlation length is approxiiftately a = 30 cm which constrains 

the frequency to f =- — » 1 kHz. At greater depths in the mixed region larger 

z fra 

values of a would relax this limitation. Brownlee (1973) feels that the require- 
ment ka » 1 is too strong. 

We have not yet considered direction of propagation. Since a is a function 
of depth, the fluctuations will vary with angle of inclination of propagation. In 
principle the solution to the problem can be approximated by integration 

q = K r 0 (z) a (z) dR 

A z '.'0 H 



as in a ray diagram calculation. In practice, although we have equation (3.40) 
for a(z) 



a(z) 



= 0 . 2 e 



z/40 



(3.40a) 



our frustration at generally not knowing how 0 varies with time, or location is 
now compounded by the non-trivial requiremeitt that 0 be known at all depths 
along the ray path. 

We have also assumed that C ( 4 ) is constant during the time of the experi- 
ment. Any inconstancy will, of C(Mjrse, cause a variation in the signal. Some 
idea of the variation in x ^ and a is possibly indicated by the scatter in data 
points in Fig. 3. 

Throughout this section, the assumption has been made, for simplicity, 
that 0 and C (^ ) are dependent only on the temperature microstructure. In 
fact, ^tudie^ of fluctuations of phase have shown that, within approximately 
5 to 10 m of the surface (depending on wind speed), fluctuations of the speed 
of sound and values of 0 are more dependent on variations of bubble radius 
and number than on temperature fluctuations. This statement is particularly 
true at sound frequencies near 60 kHz, which corresponds to the resonance 
frequency of a large population of bubbles. However, for any frequency less 
than 100 kHz, the complete description of q within the upper 10 meters of 
the ocean will require information about bubbles as well as the temperature 
structure. __ 



Finally, we must remind the reader that this report deals only with 
microstructure. To the extent that our study is correct we have evolved 
a way to understand and to predict the sound fluctuations that are universal, 
and repeatable. These fluctuations are in fact the minimum fluctuations to 
be experienced at sea. Any larger scale fluctuations due to internal waves 
or multipath propagation will be superimposed on the ones that we have 
studied so far. A fertile field of study remains for the curious. The inverse 
problem, to predict the cr from a knowledge of a or now appears also to be 
a timely subject for research. ^ 



4. Conclusion 



A truncated "universal" spatial correlation function of the excess index of 
refraction has been derived for the upper ocean. The correlation function is 
predictable from a knowledge of only the depth of the acoustic experiment. 
Although its applicability is appropriate only for reasonably well-mixed water, 
and there only for small spatial lags, it appears to be essentially correct to 
spatial copelations beyond the correlation length, that is, beyond correlation 
values e . A Gaussian extrapolation of this universal correlation function 
permits previously derived theoretical expressions for sound phase and ampli- 
tude fluctuations to be used without further changes. The predictive technique 
is simple, and appears to work well for a large range of experimental conditions. 
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